An equilibrium problem for the limiting eigenvalue 
distribution of rational Toeplitz matrices 

Steven Delvaux* Maurice Duits^ 

Abstract 

We consider the asymptotic behavior of the eigenvalues of Toeplitz matrices with rational 
symbol as the size of the matrix goes to infinity. Our main result is that the weak limit of the 
normalized eigenvalue counting measure is a particular component of the unique solution to 
a vector equilibrium problem. Moreover, we show that the other components describe the 
limiting behavior of certain generalized eigenvalues. In this way, we generalize the recent 
results of Duits and Kuijlaars [7] for banded Toeplitz matrices. 

Keywords: Toeplitz matrix, rational function, generalized eigenvalues, (vector) poten- 
tial theory. 



1 Introduction 

For an integrablc function / on the complex unit circle {z 
of size nxnis denned by 
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where fk is the kth Fourier coefficient of / 
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1} the Toeplitz matrix T n (f) 
\ 

(1-1) 

/ nxn 

(1.2) 



The function / is called the symbol of T n (f). In this paper we will be interested in symbols / 
that are rational. That is, we assume that there exist polynomials A, B\ and B 2 such that 



/(*) = 



A(z) 



B 1 (z)B 2 { Z y 



(1.3) 



where the roots of B\ (or B2) lie inside (or outside) the unit circle. Thus we do not allow / to 
have poles on the unit circle. We take A so that it has no common roots with B\ and B%. 
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Note that if B\{z) = z q , q € N, and B%(z) = 1, then (jl.3p reduces to the Laurent polynomial 



where 

p = deg A — q. 

Thus we have = for all fc > p and for all k < —q. The matrix T„(/) is then a banded 
Toeplitz matrix. The integers p and q in (|1.4|) correspond to the outermost non-zero diagonals 
in the lower and upper triangular part of this matrix, respectively. For a detailed discussion of 
banded Toeplitz matrices see [2]. 

We are interested in the asymptotic behavior of the eigenvalues of T n (f) as n — > 00. It 
is known that the eigenvalues accumulate on a particular curve in the complex plane that we 
will introduce shortly. Moreover, there exists a measure on this curve describing the limiting 
distribution of the eigenvalues. It was shown in U that for banded Toeplitz matrices this limiting 
distribution is subject to an equilibrium problem that is naturally constructed out of the symbol. 
The purpose of the present paper is to extend this result to the case of rational symbols. 

Let us first review some results on the asymptotic behavior of eigenvalues of rationally gen- 
erated Toeplitz matrices. Let / be as in Q1.3[l and T n (f) the associated Toeplitz matrix. Denote 
the spectrum of T n (f) as 

sp T n (f) - {A G C I det(T n (/) - XI) = 0}. 

To describe the asymptotic behavior of the spectrum we introduce, as in [13) . two different 
limiting sets 

liminf sp T n (f) 

n — >oo 

consisting of all A G C for which there exists a sequence {A„}„ e N, with A„ € sp T n (f) converging 
to A, and the set 

limsupsp T n (f) 

n — >oo 

consisting of all A for which there exists a sequence {A„} ng N, with A„ € sp T n (f) having a 
subsequence converging to A. 

It turns out that these limiting sets can be described in terms of solutions to the equation 

o=m-\= (1.5) 

B 1 (z)B 2 {z) 

where 

A x (z):=A(z)-\B 1 (z)B 2 {z). (1.6) 
Following the analogy with (| 1 .4|) , we define 

g:=degBi, p := max(deg J 4,degi?i degi? 2 ) - <?■ (1.7) 

To avoid trivial cases, in what follows we always assume that p, q > 1, see e.g. [B]. We also 
assume without loss of generality that 

gcd{fc I f k + 0} = 1, (1.8) 

see [1 Page 263]. 
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Note that A\(z) in (|1.6[) is a polynomial of degree p + q in z, with each of its coefficients 
depending linearly on A. There can be at most one value of A € C for which the leading coefficient 
vanishes. For all other A £ C, the polynomial A\{z) has precisely p + q roots z = z(X) (counting 
multiplicities) and we label them by absolute value as 

< ki(A)| < \z 2 (A)\ < ... < \z p+q (X)\. (1.9) 

In case where two or more subsequent roots in (jl.9l) have the same absolute value, we may 
arbitrarily label them so that (|1.9p is satisfied. For the special value of A for which the polynomial 
A\ has less than p + q roots, say p + q — k, we again order these roots 21(A), . . . , z p + q -fc(A) as in 
(|1.9p and then we set z p + q -k+i(X) = • ■ • = z p+q {X) = 00. 
Define the curve 

r :={AeC||z 9 (A)| = |z 9+1 (A)|}. (1.10) 

The fact of the matter is that 



liminf sp T n (f) = limsupsp T n (f) = T a . (1-H) 

n >oo n — toc 

This result was first established by P. Schmidt and F. Spitzer [T3] in the banded Toeplitz case 
(|1.4p . using a determinant identity by H. Widom [IB]. The generalization to the case of rational 
symbols (|1.3[) is due to K.M. Day [BJ, based on an extension [5] of Widom's determinant identity. 

Let v n be the counting measure on the eigenvalues of T n (f) 

v n =- y S x (1.12) 

Aesp T n (f) 

where S\ is the Dirac measure at A and each eigenvalue is counted according to its multiplicity. 
It turns out that the measures v n converge weakly to a measure /io on To. 

In the banded case (| 1 .4|) the measure /io is known to be absolutely continuous, and an explicit 
expression for this measure was given by I.I. Hirschman [8 . An alternative representation of /io 
can be obtained by setting k = in p,14[) below, cf. [7] . Further results about no in the banded 
case can be found in [21 [3 [3 [15] . 

For Toeplitz matrices with rational symbol (|1.3|) the limiting eigenvalue measure does not 
need to be absolutely continuous. Indeed, it was shown by Day [B] that this measure has an 
absolutely continuous part together with at most two point masses. 

Finally, we turn to the results of [7]. Consider the general system of curves 

r fc = {A e C I \z q+k (X)\ = |z g+fc+1 (A)|}, (1.13) 

for k = —q + 1, . . . ,p — 1. Each curve Tk consists of finitely many analytic arcs. We equip every 
analytic arc of with an orientation and we define the +-side (or — side) as the side on the left 
(or right) of the arc when traversing the arc according to its orientation. 
For k = —q + 1, . . . ,p — 1 we define the measure 




(1.14) 



on the curve T^. Here dX denotes the complex line element on each analytic arc of Tfc, according 
to the chosen orientation of r&. Moreover, Zj + (X) and zj_ (A) are the boundary values of Zj(X) 
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obtained from the +-side and — side respectively of These boundary values exist except for a 
finite number of points. Note that (|1.14[) is actually independent of the choice of the orientation. 

For the banded case it is shown in [JJ that each /x^ is a finite positive measure. Moreover, fj, 
is the measure of Hirschman, that is, the limit of the normalized eigenvalue counting measures 
v n as given in (|1.12p . The main observation in [7J is that the system of measures {[i k } k together 
uniquely minimizes an energy functional defined on the system of curves {T k } k . 

The purpose of this paper is to prove that also for rational symbols the measures {/J. k } k 
minimize an energy functional, thus generalizing the results in [7] . The general definition of the 
energy functional involves point sources that do not occur in the banded Toeplitz case. This is 
related to the phenomenon that the limiting eigenvalue distribution possibly has point masses 
for rationally generated Toeplitz matrices, as mentioned before. We also emphasize that the \x k 
are absolutely continuous. It is to be understood that fi$ is the absolutely continuous part of 
the limiting eigenvalue distribution, with the possible point masses removed. Our results will be 
stated in detail in the next section. 

2 Statement of results 
2.1 Auxiliary definitions 

First we introduce some definitions that will be used in the statement of our main theorems. For 
k = these definitions will be essentially the ones of Day [5], but we will state the definitions 
for general values of k G {— q + 1, . . . ,p— 1}. 

Definition 2.1. Recall the notations (|1.6|) - (|1.7[) . Define the coefficients a k ,b k E C, k = 

-q, by 

A x (z) = A(z) - XB 1 {z)B 2 (z) =: £ (a* - Xb k )z"+ k . (2.1) 

k——q 

(Note the index shift.) Define Ai, A2 G C := C U {°°} such that 

a_ g — = 0, and a p — b p X 2 = 0, 

respectively. Define k\,k 2 £ {1, . . . ,p + q} to be the maximal indices for which 

a k - b k \x = 0, k = -q, . . . , -q + ki - 1, 

and 

a k -b k \ 2 =0 1 k=p-k 2 + l,...,p, 
respectively. Here we make the convention that a k — b k Xj = if b k = and Xj = 00, j = 1, 2. 

The numbers Ai , A2 in Definition 12.11 are the unique A- values for which the polynomial Ax 
has some of its roots equal to (for Ai) or to 00 (for A2). In fact, the numbers k\ and k 2 are 
chosen such that A\ 1 has k\ roots equal to zero, and A\ 2 has k 2 roots at 00. For all other values 
of A, A\ has precisely p + q roots (counting multiplicities) which are all non-zero and finite. 

Remark 2.2. By definition we have that k\, k 2 < p + q. The case where k\ > p + q or k 2 > p + q 
cannot occur since it would imply that all coefficients of A\(z) are equal up to multiplication 
with a scalar. This would then imply that the numerator and denominator in p.3[) are equal 
up to a scalar factor, contrary to our assumptions. Note also that it is possible to have cither 
k\ = p + q or k 2 = p + q, but not simultaneously. Indeed, in the latter case we would have that 
A\(z) — (a- q — A6_ 9 ) + (a p — Xb p )z p+q and hence Bi(z)B 2 (z) = b- q + b p z p+q . The latter implies 
that either p = or q = which contradicts the assumption p, q > 1 made in the introduction. 
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Definition 2.3. For each k = 



—q + 1, .. . ,p — 1 define 



and 



q + k , 

mi,fc = max [ 1 — ,0 ] € [0, 1), 

m-2.1, = max f 1 - P fc ^ ,0 ) G [(). I.) 



m fc = 1 - mii - m 2 ,fc. 



(2.2) 
(2.3) 
(2.4) 



The numbers mii and TO 2j fc will be the weights of certain point masses, see further. The 
quantity m,k will be the total mass of the measure in (|1.14|) . Occasionally we will also consider 
mi,ki m 2,k and mfc for the indices k — — q or k — p. 

Note that the m& are strictly positive for all fc. Indeed, from the definition of rni.fc, m 2j fc, 
and the fact that k\ < p+q and fe 2 < p+<7 it is easy to check that rrik > for k = —q+1, . . . ,p— 1. 
Moreover, if = for some /c then k± = fc 2 = p + q. However, in Remark 12. 21 we observed that 
this is not possible. 

Example 2.4. Consider the banded case, i.e. bo = 1 and 6fe = for all other k. In that case we 
have Ai = A 2 = oo, k% = q and fc 2 = p. Then the numbers TO^fc, m 2! / c and mfc in Definition! 
are given in the following table 
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(2.5) 



The last row of (|2.5p contains the masses of the measures /Ltfc appearing in [7] . 
Example 2.5. Here are two examples of possible behavior when q = 4 and p = 2: 
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(2.6) 



for the case where (fci, fc 2 ) = (3, 1) or (5, 3), respectively. This occurs e.g. for the rational symbols 

. l + 2 2 +2z 6 




l+z^+z 3 +z'- + z e 
l+z 2 +z 3 +2z''+2z fi 

l + Z 2 +Z 3 +Z 5 +2 6 



(left table), 
(right table). 



2.2 The equilibrium problem 



Below we will consider measures fi supported on contours in C. If the support is unbounded 
then we will assume that 



log(l + \x\) dfi(x) < oo. 
For such a measure [i define its logarithmic energy as 



log 



1 



V - y\ 



dft(x) dfi(y). 



(2.7) 
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Similarly, for measures /i, v define their mutual energy as 

Ifav) = [[ log^^ dii{x) dv(y). (2.8) 
JJ \x-y\ 

Definition 2.6. (Compare with [7 ) We call a vector of measures v = (i/_ g _|_i, . . . , v p -\) admissi- 
ble if Vk has finite logarithmic energy, vp. is supported on T^, and Vf. has total mass Vkfik) = m k 
for every k = —q + 1, . . . ,p — 1, recall (|2.4p . 

The energy functional J is defined by 

p-l p-2 

J ( i7 )= Y 7 K)~ I{v k ,V k+ i) 
k=-q+l k=-q+l 



Here we define XAi^oo to be f when Ai ^ 00 and otherwise. The quantities xa 2 #oo, Xfei<p+g 
and Xk 2 <p+q are similarly defined. 

The equilibrium problem is to minimize the energy functional (|2.9p over all admissible vectors 
of positive measures v. 

The equilibrium problem may be understood intuitively as follows. On each of the curves 
Ffc one puts charged particles with total charge ra k . Particles that lie on the same curve repel 
each other. The particles on two consecutive curves attract each other, but with a strength that 
is only half as strong as the repulsion on a single curve. Particles on different curves that are 
non-consecutive do not interact with each other in a direct way. In addition, if Ai 7^ 00 and 
k\ < p + q then we have an external field acting on the particles on the curve T^ q+ k! ■ Similarly, 
if A2 7^ 00 and &2 < P + q we have an external field acting on the particles on r p _fe 2 . The external 
fields come from point charges at A = Ai and A = A2, respectively. The minus signs in (|2.9p imply 
that these point charges are attractive. Such external fields are sometimes referred to as 'sinks'. 
Note that there are no external fields acting on the other measures Uk, k £ {—q + k\,p — /C2}. 

Note that the external fields acting on the measures i>- q+ k 1 and v p -k 2 do not occur in [7]. 
Indeed, in that case we have Ai = A2 = 00 and hence XAi^oo and xa 2 #oo in (|2.9[) vanish, see 
Example 12.41 

For more information on equilibrium problems with external fields, see |10l 1111 IT2"] . 

Remark 2.7. In order for the above equilibrium problem to make sense, we need the energy 
functional J in (|2.9p to be bounded from below. A proof of this boundedness will be given in 
Lemma 13.71 For the boundedness it is important to note that 

Ai g r_ 9+fel and A 2 g r p _ fe2 , (2.10) 

which follows immediately from the definitions of k% and /c2- Hence the sinks are not on the 
contours on which they are acting. 

The following is our main theorem. 

Theorem 2.8. Recall the notations in (jl.7p . (jl.9p . (|1.13p and (|1.14p and assume that p,q > 1. 
Then 

(a) The vector of measures fl = (A*fe)fc = 1 L 9+ i defined in (|1.14p is admissible. 



lo S 7Z T7 dv ~q+ki (x) 



ki J |x — Ai| 

XX^ooXk 2<P+q f lQ 1 d {x y (29) 

\x - A 2 



G 



(b) For each fc G {— q + 1, . . . ,p — 1} there exists a constant fegl such that 

logfr-- — r d Vk(x)- [ lo gTT~ — i d(jL k+ i(x)- [ log 1 d^ k -i(x) 
\X-x\ J |A — x| 7 |A — re] 

h l0g ^A^ fa l0g p^=^' (2 - 11} 

for A € Tfe \ {Ai, A2}. -ffere we Zet /i^ q and /i p be the zero measures. 

(c) pi = (Mk)k=-g+i is i/ie unique solution to the equilibrium problem described above. 

Theorem 12.81 will be proved in Section [3] Note that the equalities in Part (b) are precisely 
the Euler-Lagrange variational conditions of the equilibrium problem. Part (c) will then be a 
consequence of the convexity of the energy functional J and the fact that J is bounded from 
below. 



2.3 The measures /j,k as limiting distributions of generalized eigenval- 
ues 

It was proved in [7] that in case of banded Toeplitz matrices, the measures (x^ for k 7^ also 
have an interpretation of being the limiting measures for certain generalized eigenvalues. For the 
rationally generated Toeplitz matrices such a result remains valid. 



Definition 2.9. For k = —q + 1, . . . ,p — 1 and n > 1 we define the polynomial Pk. 

P fe ,„(A) = detT„(z- fe (/(z)-A)) 
and we define the kth generalized spectrum of T n (f) by 

sp fc T„(/) = {A e C I P fe , n (A) = 0}. 
Finally, we define [ik,n as the normalized zero counting measure of sp fc T„(/) 

Mfe.n = — Sx, 
n 

Aes Pfc T„(/) 

where in the sum each A is counted according to its multiplicity as a zero of Pk, n - 



by 



The Toeplitz matrix T n (z k (f — A)) in Definition 12.91 may be interpreted as a shifted version 
of T n (f - A) = T n (f) - XI. For example if k = 2 then by (TOD we have 



T„(z- 2 (/-A)) 
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(2.12) 



We will show that for each k the sequence {fik,n}n has a limit. Moreover, the limiting measure 
will have point masses at A = Ai (if Ai =/= 00) and A = A2 (if A2 ^ 00) with weights at least m^t 
and TO2,fc, respectively. On the other hand, if Ai = 00 or A2 = 00 then the total mass of the 
limiting measure is reduced with at least mx t k or m2,k respectively. These facts can already be 
seen at the level of the finite-n measures fj,k,n as the following proposition shows. 
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Proposition 2.10. Let k G {— a + l,...,p— 1}. Then the polynomial Pk, n {X) satisfies the 
following properties. 



(a) 



(b) 



P k ,„(X) is divisible by (A - Ai)™ 1 -*™^. if Xi G C, 

■Pfe,n(A) has degree at most (1 — m^jn + c, if Ai = oo, 

P fc , n (A) is divisible 6y (A - A 2 ) m2 - fc ™" c , z/ A 2 G C, 

-Pfe,n(A) /ias degree at most (1 — m 2j fc)n + c, i/ A 2 = oo, 



where c > is a constant depending only on the symbol f . 

Denote by Qk.n{X) the quotient polynomial obtained from -Pfc, n (A) by removing all its factors 
(A — Ai) (if Ax S Cj and (A — A 2 ) (z/ A 2 G CJ. TTien we have that 

(c) Qk,n{X) has degree at most m^n + 2c. 

Proposition 12. 101 will be proved in Section B~T1 

Since the measure (|1.14j) is absolutely continuous, the best one can hope for is pk being the 
absolutely continuous part of the limiting fcth generalized eigenvalues distribution. This means 
that the possible point masses at Ai and A 2 should be stripped out in /i^. This turns out to be 
indeed the case. 

Theorem 2.11. Let k G {-q + 1, . . . ,p - 1}. Then 

liminf sp fc T„(/) = limsupsp fc T n (/) = T fe , 

n >oc n — > qq 

and 

lim / 0(A) dfik.n(X) = / </>(A) d/x fc (A) + XA^ooTOi.fe^CAi) + XA 2 ^oo™2.fc0(A 2 ) (2.13) 
n ^°°Jc Jc 

holds for every bounded continuous function (f> on C. 

From (|2.13j) we see that and to 2j /c are the weights of the point masses at Ai and A 2 in 
the limiting fcth generalized eigenvalues distribution, if present. 



2.4 Organization of the rest of the paper 

The rest of this paper is organized as follows. Section contains the proof of Theorem 12.81 In 
Section [H we prove Proposition 12.101 and Theorem 12.111 Most of the proofs are inspired by the 
proofs given in [7] for the corresponding statements in the banded case, hence we will often refer 
to that paper. Finally, some illustrations of our results are given in Section [5l 



3 Proof of Theorem 

3.1 Proof of Theorem [23](a)-(b) 

In this section we will prove Theorem !2.8f aWb). First we recall some elementary definitions 
and properties involving the algebraic equation A\(z) = 0. 

Definition 3.1. ([U Section 11. 2], [7]) A point A G C is called a branch point if A\(z) = has 
a multiple root. A point A G Tk is an exceptional point of Tk if A is a branch point, or if there 
is no open neighborhood U of A such that Tk PI U is an analytic arc starting and terminating on 
dU. 
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Proposition 3.2. Let k G {—q + 1, . . . ,p — 1}. Then the set T k in (|1.13p is the disjoint union 
of a finite number of open analytic arcs and a finite number of exceptional points. The set T k 
has no isolated points. 



The proof of this proposition is similar as in 2, Theorem 11.9], ff\. [13{. The condition (|1.8|) 
is needed to ensure that the Tk are proper curves, i.e., they are 1-dimcnsional subsets of C. 
A major role is played by the functions w k which, for k = —q + 1, . . . ,p — 1 are defined by 

q+k 

w k (\) = Y[z 3 (X), XeC\r k . (3.1) 

The function w k is analytic in C \ T k . Occasionally we will also consider w k for the indices 
k = —q or k = p. 

Note that (|1.14[) may be written alternatively as 

dfJ-k(X) = ~ '■ 7TT 7TT d\. 

2m \w k+ (X) Wk-(X)J 

To discuss the integrability of this measure, we will need the asymptotic behavior of ^ ^ . The 
relevant facts are listed in the following proposition. 

Proposition 3.3. Let k G {— q+ 1, .. .,p — 1} and recall the notations in Section \2.1\ Then the 
following statements hold 

(a) For any Xq G C \ {Ai, A2}, there exists an I G N := {1, 2, 3, . . .} such that 

<(A) f 0((A - A o )- 1+1 /0, if A G C 
u> fc (A) 1 ©(A- 1 " 1 /'), «/A = co 

as A — > Ao wtt A G C \ Tfc . We have I = 1 unless Xq is a branch point. 

(b) Assume X\ 7^ A2. Then near the point X\ there exists an I G N such that 

w' k (X) j 1 - T ^f + 0((X-X 1 )-^/ 1 ), ifX 1 eC 

1 — mi 



OO 



as X — > Ai im£/i A G C \ Tfc. 

(^cj Assume X\ ^Xi. Then near the point A2 there exists an I G N suc/i i/iai 

w' k (X) = f ^ + O((A-A 2 )- 1+1 /0, z/A 2 GC 

^fe(A) \ !^ + (A- 1 - 1 /') ; z/A 2 = oo 

as A — > A2 tuit/i A G C \ T k . 
(d) Lf Ai = A2. then 

^fc(A) \ ^ + 0(A- 1 - 1 / ; ), z/A 1= oo 
as A — > Ai = A2 with X G C \ T k . 
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Proof. First we make some general observations. For any A G C \ {Ai,A2} the polynomial 
A\(z) = has roots Zj(X), j — 1, . . . ,p + q, all of which are finite and non-zero (although some 
of the roots might occur with higher multiplicity) . Zero roots or infinite roots can only occur if 
A G {Ai, A2}. As A — > Ai, then the k\ smallest roots z,(A) tend to zero like a power (A — Ai) 1 '* 51 , 
while the other roots converge to non-zero constants. Similarly, as A — > A2, then the ki largest 
roots Zj(A) tend to infinity like a power (A — X^ 1 ^ 12 , while the other roots converge to non-zero 
constants. These facts all follow from Definition 12.11 

To prove the first equality of Part (a), fix A G C \ {Ai, A 2 } and denote with Zj = Zj(Xo), 
j = 1, . . .,p + q, the roots of A\ {z); by the discussion in the previous paragraph we have 
Zj € C \ {0} for each j. Pick one of the roots Zj which has the highest multiplicity I. Writing 

A x (z) = (z- Z] ) l P(z) + (A - Xo)Q(z) 
where P(z) and Q(z) are polynomials with P(zj) ^ and Q(zj) 7^ 0, it follows that 

Zj(X) = Zj +c(X- X Q ) 1/l , A^A , 
for some constant c. By taking the logarithmic derivative we obtain 

fM^A-Ao)-^), A-Ao. 

Zj{A) 

The first equality in Part (a) follows from this and the fact that 

<W _y Z 'M) 

The second equality in Part (a) (for the case Ao = 00) is proved in a similar way, this time 
using a decomposition 

A x {z) = {z - Zj) l XP(z) + Q{z) 

where again P(zj) ^ and Q(zj) 7^ 0. 

To prove Part (b), first assume that Ai 00. From the discussion in the first paragraph of 
this proof we obtain 

^ iAi ' +0((A-A 1 r 1+1 /<-), A -Ax, 



Z j(X) fei(A-Ai) 
for all j = 1, . . . , kt, while 

z'(X) 



0((A-A!)- 1+1 /'), A — > Ai, 



for j = ki + 1, . . . ,p + q and a suitable I S N. Hence 

'q + k \ 1 



k\ J X — Ai 
1 - mi, fc 
A-Ai 



+ 0((A-A!)- 1+1 / F ) 



+ 0((A-A 1 )- 1+1 / / ) 



by virtue of (|2.2p . with / = max(/,fci). Similarly one can prove the case where Ai = 00. The 
proofs of Parts (c) and (d) are similar as well. □ 
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Proposition 3.4. For each k we have that fik in (1 1 . 14f> is a positive measure on Tk with total 
mass /ife(r^) = rrik as defined in (|2.4p . 

Proof. First we prove that the density (|1 . 14[) is locally integrable around the points Ai, A2, 00 
(at least those of them which lie on the curve Tfc). For 00 this follows from the second equality 
in Proposition 13. 3f a). For Ai this follows from Proposition 13 . 3f b) and taking into account that 
the 1/(A — Ai) terms at the +-side and — side in (|1.14p cancel; a similar argument holds for the 
point A2. 

The fact that the measure /!& is positive follows from a Cauchy-Riemann argument as in [7j 
Proof of Proposition 4.1]. 

Finally, the statement that fj,k(^k) = f^k follows from a contour deformation argument as in 
[7J Proof of Proposition 4.1]. More precisely, we have 

Mfc( r fe) : = 7T- 77V 7TT dX 

2ni J Fk \w k+ {\) w k -{X)J 

1 f wUX) ^ n fw'JX) \ n (w'JX) , , \ 

2iriJ c Wk(X) \Wk(X) ) r \Wk{X) J 

(3.2) 

where C is a clockwise oriented contour surrounding Tk and those points Ai and A2 which are 
finite, and where Res(/i, A) denotes the residue of h at A. Note that <\'6.'2\i is valid even when some 
of the Xj lie on the curve r^, j = 1, 2, thanks to the local integrability of Hk around these points. 
Applying the residue theorem once again, this time for the exterior domain of C, we then find 
for the first term in (|3.2p that 

1 f<W „ D „/X(A) 



, . dX = —Res — ^tt, A = 00 ] . (3.3) 

2mJ c w k {X) \w k (X)' J y ' 

The fact that /J,k{^k) = m,k then follows from (|3.2p - (|3.3p and the residue expressions in Propo- 
sition [331 recall also (|2~4]l . □ 

Proposition 3.5. For each k we have that 

dfj, k (x) = — ^— + xa 1 #oo-t — : — hxA 2 #ooT — i/AeC\r t (3.4) 



A — x Wk(X) X — Ai A — A2 



log I A - x\ dfik(x) = -log|i«fc(A)| + XA 17 too(l - mi,fe)log|A - Ai| 

- XA 2 ^ooTO2,fclog|A - A 2 | +a kl ifXeC (3.5) 

for a suitable constant ctk ■ 

Remark 3.6. The quantities in the above proposition are only well-defined if A 7^ Xj, j = 1, 2. 
However, one easily checks that Ai and A2 are removable singularities for the right hand sides of 
both (|3.4p and (|3.5p . due to the continuity of the corresponding left hand sides. 

Proof of Proposition \3.5l The proof of (|3.4p follows by contour deformation in a similar way as 
in the proof of Proposition ^. 41 The relevant expression is now 

J Tk X - x 2m J Tk X - x \Wk+{x) W k -[x)J 2m J e X — % Wk(x) 



, Xa^oot r-Res '- r ,x = Ai +xa 2 #oot r-Res — -. r, X = X 

Wk(X) r X-Xi \Wk{x) J A-A 2 \w k {x) 
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where C is a clockwise oriented contour surrounding Y k , the point A, and those points Ai and 
A2 which are finite. Now the integrand in the integral over C has zero residue at infinity and 
therefore this integral vanishes. Using the residue expressions in Proposition ^ . 31 one then arrives 
at the right hand side of (|3.4p . Finally, the proof of (|3.5|) then follows by integrating (|3.4p . see 
also 0. □ 

Now we are ready to prove Theorem I2.8f a)-(b) . 

Proof of Theorem \2.8\( a). Taking into account Proposition 13.41 it suffices to show that the loga- 
rithmic energy I(fJ, k ) is bounded for each k 6 {—q+1, . . . ,p— 1}. The latter follows by integrating 
(|3.5[) over /ifc(A). Then the left hand side becomes — I(p k ), so it suffices to show that each of 
the four terms in the right hand side is bounded. For the fourth term this is evident since ix k 
has finite mass. For the two middle terms this follows from our earlier observation that jj, k is 
integrable around Ai and A2 (assuming they are on the curve T^), which is still true when mul- 
tiplying with the logarithmic singularities log |A — Ai | and log |A — A2 1. A similar argument holds 
for the first term. □ 

Proof of Theorern \ 2. t^ b). The proof of Part (b) follows from (|3.5[) and the auxiliary results 



-log|w fe+ i(A)|+21og|w fe (A)|-log|w fe _i(A)| = log 
for A e T k \ {Ai, A 2 }, and 



nL^-(A) 2 



n:^(A)n^(A) 



= log 



Zfe+i(A) 



= 



J -1/h k=-q+ki , , 

-m l , k+l + 2m x , k -m 1 , k . l = ^ Q fc€{ _ ?+1><i . jp _ 1}u _ 9 + fel} , (3-6) 

/ -l/k 2 k = p- k 2 

-m 2 , fc+ i+2m 2)fe -m 2)fc _i = j Q g {— g + 1, . . . ,p — l}\{p — k 2 } (37) 

Here the boundary terms mi &, m 2j t for = —5 or fc = p are defined by the usual formulae 
(|2.2p - (|2.3[) . These considerations imply the desired result for A € \ {Ai, A 2 }; the cases A = Ai 
and A = A 2 then follow by continuity. □ 

3.2 Proof of Theorem &Mc) 

To prove Theorem 12. 8f c) we rewrite (|2.9p in the following way, compare with Eq. (2.12)]: 

■ / " 7 )=l L '[ — - — ) I + -r^—I{v- q+kl ) + — ^/(^_ fc2 ) 

1 — ' ' 1 1 2k x m kl 2k 2 m k2 

dv- q+kl (x) - XA ^°° y log r; J 1 ^! dv p _ k2 {x). (3.8) 

We leave it to the reader to check the correctness of this identity; note that the calculation makes 
use of the auxiliary result 




-mfe+i + 2m k - m k ^ x 



1/fci k = -q + k 1 ^p—k 2 

l/fca k=p — k%^—q + k\ 

l/fci + l/fc 2 k = —q + k 1 =p~k 2 

fc S {-g+l,...,p- l}\{-g + fci,p-fc 2 } 
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for k e {—q+1, . . . ,p— 1}, which follows from (|2.4p and (|3.6p (|3.7p . Here we recall the boundary 
values m_ g = m p = 0. 

We also invoke the fact that 

- v 2 ) > 0, (3.9) 

whenever v± and v 2 are positive measures with v\(C) = ^(C) < oo. This is a well-known result 
if v\ and v 2 have bounded support [T^] . If the support is unbounded this is a recent result of 
Simeonov [Ti] , 

Lemma 3.7. The energy junctional (|2.9p is hounded from below on the set of admissible vectors 
of measures v. 

Proof. From (|3.8p ~ (|3.9p we see that in order to show that the energy functional J{v) is bounded 
from below, it is sufficient to show that 

— L-I(„_ q+kl ) - 2^ f log^— dv- q+kl {x) (3.10) 
2kim kl k x J \x-\i\ 

and 

— — I(f p -k 2 ) ~ A ff°° / log | r-r dv p - k2 (x) 

2k 2 m k2 k 2 J \x - \ 2 \ 

are both bounded from below on the set of admissible vectors of measures v. Let us check this 
for the first term (|3.10p . We will use that 



q+ki 



a fact already observed in (|2.10p , which follows immediately from the definition of ki . Now we 
distinguish between two cases. The first case is when Ai = oo. Then the second term in (|3.10p 
drops out while on the other hand oo = Ai T- q+kl , so the contour T- q+kl is bounded and 
therefore the first term in (|3.10p is bounded from below as well. 

The second case is when Ai ^ oo. Then standard arguments from potential theory show 
that the expression (|3.10p is minimized precisely when v- q+kl is the balayage of the Dirac point 
mass at Ai onto the curve r - q + kl , and in particular this expression is bounded from below [121 
Chapter 2]. □ 

Remark 3.8. The above proof goes through because the constant factor in front of the first 
term in (|3.10p is precisely 1/2. If this constant factor is different from 1/2 then the connection 
with balayage measures breaks down, and in fact if the constant is larger than 1/2 then the 
energy functional is not bounded from below anymore. 

Proof of Theorem ] 2.81( c). Assume that pi is a vector of admissible measures satisfying the equal- 
ities in Theorem 12. 8f b). and let v be any admissible vector of measures. We need to prove that 
J{y) > J(/2) with equality if and only if v = fx. Note that the equalities in Theorem 12. 8T bl are 
precisely the Euler-Lagrange variational conditions of the equilibrium problem. The result then 
follows from the fact that the energy functional J is convex and bounded from below. More 
precisely, one can use exactly the same argument as in [71 Proofs of Lemma 2.3 and Theorem 
2.3(c)], taking into account (|3.8p - (|3.9p . There are some modifications induced by the external 
fields, but this does not lead to problems since the latter act in a linear way on the measures. □ 
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4 Proofs of Proposition 12.101 and Theorem 12.111 



4.1 Proof of Proposition 1^.101 

The proof of Proposition ^. 101 is based on the reduction of a rationally generated Toeplitz matrix 
into banded form, which will then allow us to follow the proof in [JJ Proof of Prop. 2.5]. Let us 
recall from (|1.5[) that 

f(z)-X= ^) 
M ' B 1 {z)B 2 {zY 

where the numerator A\{z) is a polynomial in z. Then we claim that for any k G {— q + 
l,...,p — 1} and for any n sufficiently large, the rationally generated Toeplitz matrix with 
symbol z~ k (f(z) — A) can be reduced into banded form by the factorization 

L n T n (z- k (f(z)-\))Rn=T n (z-«- k A x (z))+ ( C fj , (4.1) 

V / nxn 

where 

L n = T n (B 2 (z)), R n := T n (z^B^z)) , 

are non-singular lower and upper triangular Toeplitz matrices respectively. The middle factor 
in the left hand side of (|4.ip is our rationally generated Toeplitz matrix of interest, and (|4.1|) 
shows that it can be reduced to the banded matrix pencil in the right hand side. Here C is a 
matrix whose size and entries are independent of n but depend only on the symbol f(z). For 
more information on factorizations of the type (|4.ip see e.g. [H Prop. 2.12] and also [5j[9]. 
From 1|4.1[) we immediately deduce that 

P k , n (X) := detT n (z- k (f(z)-\)) 

I n n\ \ 

(4.2) 



where k is a numerical constant, given by the product of the diagonal entries of the two 
triangular factors L n and R n in (|4.1|) . 

We are now ready for the proof of Proposition 12.101 The proof will follow by expanding the 
determinant in (14. 2|) by a basic combinatorial argument, see also [7J Proof of Prop. 2.5]. 

Proof of Proposition ^. 10\f a) . The proposition is obvious if m^jj = 0. So we will assume below 
that m^t > 0, or equivalcntly 

q + k<h. (4.3) 
First we consider the case where Ai ^ oo. By expanding the determinant in (14. 2(1 we find 

1 " 

Pk.nW = - II ( a J-^0')+fc _ A& j-T(j) + fc +Xj<|C|X7r(i)<|C] S>0')) • 

7r£S„ j=i 

Here S* n denotes the set of all permutations of {1, ... ,n}, and we denote with |C| the maximum 
of the row and column sizes of the matrix C; note that this number is independent of n. By the 
band structure it follows that we only have non-zero contributions for the permutations it that 
satisfy 

-q-k<j- ir(j) <p-k for all j = \C\ + 1, . . . , n. (4.4) 
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Denote, for tt € S n , 

= {j | j - n(j) £{-q-k,...,-q-k + k 1 - 1}}. (4.5) 

The set contains all indices j for which the (j, ir(j)) entry lies in the union of the fci topmost 
bands of the banded matrix in (|4.2[) . By assumption (|4.3|) these bands include the main diagonal 
j — 7t(j) — and by definition of ki we have that the entries in these bands are all divisible by 
(A -A X ). 

Denote the number of elements of N n in (|4.5[) by \N„\. Then obviously 

P fc ,„ is divisible by (A - \ 1 ) min * res " (4.6) 

where we minimize over all permutations tt 6 5„ satisfying (|4.4p . 

Let tt e 5„ satisfy (|4.4|) . We give a lower bound for |AT ff |. Since X}j=i(j — n U)) = we 
obtain 

n n 

E(j-tC7))+ = E(t(j)-j)+, (4-7) 
3=1 3=1 

where (•)+ is defined as (a)+ = max(0, a) for a G M. From the above definitions we also have 
that 

j - 7r(j) >-q-k, if j e , . 

i-7r(7)>-«-fc + fci, ifje{|C| + l,...,n}\JV w . 

By combining (|4.7p and (|4.8p we find that 

n n 

Hz - fc + *i)(n - |iV.|) + c < £(j - 7r(i)) + = £>(j) - j)+ < + - c. 

3=1 3=1 

Here c > is a correction term which is due to the presence of the matrix C in the top left 
matrix corner in (14. 2|) : the number c is clearly bounded from above. We then obtain 



k\ — q — k ( q + fc\ 
> : n - c = 1 — n - c = 7ni ifc n - c, (4.9) 



fci \ fci 

where we used (|2.2j) and (|4.3p . and where we put c := 2c/k\. The first statement in Proposi- 
tion [2TjJa) now follows from (|4T5j) and (g^J) . 

The proof of the second statement in Proposition 12. lOf a) (for Ai = oo) is similar to the one 
above. Now one uses that all the entries aj—^u^+k ~ ^bj—irfjj+k m the bands indexed by j G N n 
have their A-coefficient bj_K(j)+k = 0, which then yields in a similar way to (|4.6[) and (|4. 9[) that 

degP fc ,„ < n - min |AT ff | < (1 - m hk )n + c, (4.10) 

as desired. □ 

Proof of Provosition [KTW b). Similar to Part (a). □ 

Proof of Provosition \2.1(]\f c). Part (c) follows immediately from Parts (a) and (b), together with 
(|2.4p . in case where Ai ^ A2. The case where Ai = A2 can be obtained as well, by observing that 
at least one of the numbers ui\^ and m-i^ must be zero in that case. The latter follows since 
otherwise the numerator and denominator in (|1.3p are equal up to multiplication with a scalar, 
contrary to our assumptions. □ 
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4.2 Proof of Theorem I2~TTT1 

To prove Theorem 12. Ill we need to manipulate the polynomial Pk, n {X). To this end we will use 
a determinant identity by K.M. Day which we state next. 

To state the identity, we need some notations. Denote with ft and 7* the zeros of Bi(z) and 
Bi{z), respectively. Recall the notation z% = Zi{\) for the roots of A\(z). Thus 

p+q 

A\{z) = c\{{z-z t {\)) (4.11) 

<=i 
<? 

Bi(z) = c x TJ(2-A) (4.12) 

i=l 

dcg B 2 (z) 

B 2 (z) = c 2 H («-7<)» ( 4 -!3) 

where c, ci, C2 are non-zero constants. 

The following theorem was proved under some additional hypotheses by K.M. Day [5]. Other 
proofs are in [3J [9], the former of them stated under the weakest assumptions. We state the 
theorem in the form that is most convenient for our purposes. 

Theorem 4.1. (Day's determinant identity). Let k G {— q+ 1, . . . ,p— 1} and let A G C\ {Ai, A2} 
6e swc/i that all roots of A\(z) are distinct. Then 

P k , n (X) = detT„(z- fc (/ - A)) = Cs(\)(w s (\)T, (4.14) 

s 

where the sum is over all subsets S C {1,2, . . . ,p + q} of cardinality \S\ = q + k and for each 
such S we have 

w s (X) := (-l) q+k (a- q - 6_ g A) fll^wj ( 4J5 ) 
and (with S := {1, 2, . . . ,p + q} \ S) 

( Zj (X) - /3 r )( 7t - Zi(X)) 



c5(A) : =n^( A ) fe n 



(2 i (A)-« j (A))( 7 t-A-)' 
je5 1 e S,r £ R, 

jeS,teT 

with R = {1, . . . , g} and T := {1, . . . , deg ^(z)}- 

Incidentally, observe that (|4.15[) can be written alternatively as 

w s (X) = (~l) p - k (a P - b p X) Y[ Zj (X). (4.16) 

We note that in case where k G {0, . . . ,p— 1}, our formulation of Theorem l4.ll follows directly 
from the one of [5]; for the case k G {— g + 1, . . . , —1} it can be obtained from the result of [S] 
by working with the transposed matrix. 

From (|4.14j) - (|4.15[) we see that for large n, the main contribution in (|4.14p comes from those 
subsets S for which |ws(A)| is the largest possible. For A G C \ there is a unique such S, 
namely 

S = S k :={l,2,...,g + fc}. 
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Now we are ready to show that the asymptotic distribution of the feth generalized eigenvalues 
of T n (f) is described by the measure [i k , together with possible point masses at Ai and First 
we prove this at the level of the Cauchy transforms. 

Proposition 4.2. Let k € {— q + 1, . . . ,p — 1}. Then 

,. f dfj, kin (x) f dfi k (x) mi lfc m 2 , k ,. 

hm / — = / +x Ai #oot t-+Xa 2 #oot r- (4.17) 

n->oc J c A — x Jc A-x A — Ai A — A2 

uniformly on compact subsets of C\T k - 

Remark 4.3. The above proposition implicitly assumes that A ^ Aj, j = 1, 2. However one 
checks that if Aj £ C \ then Aj is a removable singularity for the right hand side of (|4.17p . 
due to the continuity of the left hand side, and then the uniform convergence still applies. 



Proof of Proposition \4-S\ As already mentioned, for large n the dominant term in Day's deter- 
minant identity Theorem 14. II is obtained by taking S — Sk '■= {1, 2, . . . , q + k}. Then we find in 
the same way as in [71 Proof of Corollary 5.3] that 

. d/J, k ,n(x) .. 1 V"^ 1 r 1 P k,nW 

lim / — — ■ = lim — > — = lim — — — — - 

A — x n— too n ^ — ' A — A, n^oo n Ph n (A) 

_<(A)_ <(A) 



, \ , , \"\ +XA!#oot r- (4.18) 

ws k W Wk[X) r A-Ai 

uniformly on compact subsets of C \ Y k , where the last transition of (14.18[) follows from (|3.1|) 
and (|4.15p . Now from Proposition 13 . 51 we see that the right hand side of (|4.18[) equals the right 
hand side of (|4. 1 7|) . The proposition is proved. □ 

Now we are ready for the 

Proof of Theorem \2.11i From the convergence of the Cauchy transforms in Proposition 14.21 we 
deduce that 

in the weak-star sense, which means that (|2.13p holds for every continuous cj> that vanishes 
at infinity. Now a priori, it is not immediate that (|2.13p holds for all bounded continuous 
functions since it is possible that /ifc jrl has mass leaking to infinity as n — > 00. However, from 
Proposition 12. 101 it follows that this cannot happen, i.e., the measures {/ifc,n}n are tight. Thus 
(|2.13p holds indeed for all bounded continuous functions. For more details see [7J Proof of 
Theorem 2.6]. □ 

5 Example 

Consider the rationally generated Toeplitz matrix with symbol 

= 2z 2 -5z + 2 = (2z-l)(z-2) (5 ' 1} 

defined on the complex unit circle. We may compute the Fourier series of this symbol explicitly 
and find 

1 1 1 1 z z 2 z 3 

* ^ ~ ' ' ' ~ 12^3 ~fe2~3i~6~12 _ 24 _ 48"" 
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So the rationally generated Toeplitz matrix T n (f) looks like 



TnU) = - 



/1/6 


1/3 


1/6 


1/12 


1/12 


1/6 


1/3 


1/6 


1/24 


1/12 


1/6 


1/3 


1/48 


1/24 


1/12 


1/6 



Equation (|1.6p now becomes 

A x (z) = 1- A(2z 2 -5z + 2), 
and (|1.7p leads to p = q = 1. The roots of A\(z) are given by 

zi, 2 (A) = ^(5A±v / 9A2 + 8A), 

and they should be labeled in such a way that |zi(A)| < (22(A)! for all A. The roots Zi(A) and 
22(A) are coalescing precisely when 9A 2 + 8A = 0, so the branch points are A = and A = —8/9. 

Since p = q = 1, there is only one relevant index k in (|1.13[) . namely k — 0. The corresponding 
set Tq is simply the line segment connecting the branch points A = and A = —8/9: 

r = {A e C I \ Zl (X)\ = |z 2 (A)|} = [-8/9,0]. 

This may be checked from a straightforward calculation. 

Definitions 12.11 and 12.31 now specialize as follows: Ai = 1/2, A2 = 0, k\ = 1, fc 2 = 2, and 
TO1.0 = 0, 7712.0 = 1/2 and mo = 1/2. Thus the limiting eigenvalue distribution of the matrix 
T„(/) for n — > 00 consists of an absolutely continuous part liq with total mass 1/2, supported 
on To = [—8/9, 0], and a singular part which is a point mass of mass 1/2 at A = 0. 

The energy functional (|2.9[) now specializes to 



I(u ) - / log 



1 



1/2| 



du (x). 



(5.2) 



So liq is the minimizer of (|5.2[) over all measures v§ on Tq = [—8/9,0] with total mass 1/2. The 
second term in (|5.2p can be interpreted as an attraction of (1q towards the point A = 1/2. 

The measure liq is absolutely continuous with density given by (|1.14p (with k = and 
p = q = 1). The density can be explicitly computed, but we will omit the result since it does not 
lead to considerable insight. We only mention that the density blows up like an inverse square 
root near both endpoints A = and A = —8/9. More precisely, it behaves approximately like 
0.28/y/\X\ near A = and like 0.10/^A + 8/9 near A = -8/9. 

Figure [T] contains a plot of the limiting density. The figure shows that there is more mass 
near than near —8/9, which is due to the attraction towards A = 1/2 in (|5.2p . 

Figure [2] shows the result of a numerical computation of the eigenvalues of T n (f) with n = 60. 
Note that approximately half of the eigenvalues is located at zero, according to Proposition ^. 101 
in fact we have c = in this case. 



We may consider the following modification of (|5.ip . 

, . l + ez 1 

f( z ) 



ez 



2z 2 -5z + 2 (2z-l)(z-2)' 



(5.3) 



where e > is some small number. It is still true that Ai — 1/2 and A2 = for any e, but for e 
non-zero we now have k\ — k^ — 1, mio = "12.0 = an d ma = 1. Thus the limiting eigenvalue 
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lambda 



Figure 1: Density of the measure hq on Tq — [—8/9,0] for the symbol (|5.ip . The density blows 
up like a square root near both endpoints —8/9 and 0. There is more mass near due to the 
attraction towards A = 1/2. 

distribution of T n (f) is absolutely continuous (without point mass), it has total mass 1, and it 
is supported on the interval Tq joining the two branch points 

-4 - 5e± 2V4 + lOe + 4e 2 

9 ' (5 ' 4) 

From the above discussions, we see that the limiting eigenvalue distribution of T n (f) is 
absolutely continuous if e > and has a point mass at the origin if e = 0. To understand this, 
note that for e > the energy functional (|2.9p contains attracting point charges at both Ai = 1/2 
and A2 = (since k\ = ki2 = 1). In the limit e — > 0, the rightmost endpoint of Tq in (|5.4[) moves 
towards the point source at A2 = 0. This causes an increasing accumulation of mass near this 
endpoint which in the limit for e = gives birth to the point mass. 
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